Spin-3/2 physics of semiconductor hole nanowires: Valence-band mixing and tunable interplay 
between bulk-material and orbital bound-state spin splittings 
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We present a detailed theoretical study of the electronic spectrum and Zeeman splitting in hole quantum wires. 
The spin-3/2 character of the topmost bulk- valence-band states results in a strong variation of subband-edge g 
factors between different subbands. We elucidate the interplay between quantum confinement and heavy-hole - 
light-hole mixing and identify a certain robustness displayed by low-lying hole-wire subband edges with respect 
to changes in the shape or strength of the wire potential. The ability to address individual subband edges in, 
e.g., transport or optical experiments enables the study of holes states with nonstandard spin polarization, which 
do not exist in spin-1/2 systems. Changing the aspect ratio of hole wires with rectangular cross-section turns 
out to strongly affect the g factor of subband edges, providing an opportunity for versatile in-situ tuning of 
hole-spin properties with possible application in spintronics. The relative importance of cubic crystal symmetry 
is discussed, as well as the spin splitting away from zone-center subband edges. 



I. INTRODUCTION 

Self-assembled semiconductor nanowires have attracted a 
lot of attention recently due to their well-defined crystalline 
structure, unique electrical and optical properties, as well 
as a promising outlook for their use as building blocks in 
nanoelectronics, nanospintronics and nanobiotechnology. 132 
Rapid developments in material science and technology have 
enabled the realization of nanowires from many materi- 
als and materials combinations, such as III-V semiconduc- 
tors (e.g., GaAs, 1 ' 22,27 M 5 , 2 ' 3 ' 10 ' 12 i nAs i,7,8,29^ H-VI semi- 
conductors (e.g., CdSe, 512 CdS, 12 CdTe, 4 - 25 ZnTe 20 ' 32 ), as 
well as Si 14 ' 24 ' 28 and Ge. 23 ' 30 ' 31 By now it is possible to 
dope nanowires 22 26 ' 29 and perform heterostructure engineer- 
ing with atomic precision. 7 ' 8,1217 This enabled the realization 
of p-type 22 ' 26 or ambipolar 29 transport characteristics, as well 
as the fabrication of nanostructures (quantum dots, 8 rods, 20 
and superlattices 9 ). Based on these material-science devel- 
opments, devices such as nanowire diodes, 13 superlattices 
for nanoscale photonics and electronics, 9 lasers, 6 ' 16 resonant 
tunneling diodes, 8 and light-emitting diodes, 1312 have been 
demonstrated. 

Another area of significant research that has emerged in re- 
cent years focuses on the study of the spin degree of free- 
dom. This is relevant for the fundamental understanding of 
spin phenomena in solid-state systems, but also for potential 
applications that utilize the spin, rather than the charge de- 
gree of freedom, in electronics. 33 34 Spin injection, manipula- 
tion and detection has been achieved with the use of magnetic 
fields and diluted magnetic semiconductors. 35 ' 36 A comple- 
mentary direction of current spin-electronics research is based 
on the quantum-mechanical coupling between spin and or- 
bital degrees of freedom. For example, the tunability of struc- 
tural inversion-asymmetry induced spin-orbit (SO) coupling 37 
could be used to manipulate spin-polarized currents in ballis- 
tic mesoscopic channels. 38 



To study SO effects in semiconductors, charge carriers from 
the valence band (i.e., holes) are particularly interesting since 
they are subject to an inherently strong SO coupling even in 
the bulk. In addition, as valence-band states are predomi- 
nantly p-like (i.e., have orbital angular momentum L = 1), 
they are characterized by a total angular momentum of either 
J = 3/2 or J = 1/2. The topmost valence bands are de- 
scribed by the former; thus they behave like spin-3/2 particles 
and, when confined in nanostructures, exhibit quite counter- 
intuitive quantum effects. For example, the spin-3/2 nature 
of the topmost valence bands leads to large anisotropies of 
the hole g factor in quantum wells, 39,40 wires, 27 point con- 
tacts, 41 ' 42 quantum dots, 43 ^ 15 and localized acceptor states. 46 

The bulk dispersions for the conduction and valence bands 
of a typical semiconductor in zero magnetic field are shown 
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FIG. 1: (Color online) Schematic dispersion of bulk-semiconductor 
conduction and valence bands (solid lines). The dashed lines depict a 
situation where the J = 3/2 valence bands are split due to quantum 
confinement, e.g., in a two-dimensional hole gas. 
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as the solid curves in Fig. 1. The four-fold degeneracy of 
the J — 3/2 valence-band edges is lifted at finite wave vec- 
tor k, 47 giving rise to separated heavy-hole (HH) and light- 
hole (LH) branches whose states are distinguished by the z- 
axis spin-3/2 projection quantum numbers J z = ±3/2 and 
J z = ±1/2, respectively. Size quantization, e.g., due to 
a quantum-well confinement, causes an energy splitting be- 
tween the quasi-twodimensional (2D) HH and 2D LH band 
edges. 39 ' 48 49 This HH-LH splitting is illustrated by the dashed 
lines in Fig. 1 . At finite wave vector k ^ for motion perpendic- 
ular to the quantum-well growth direction, a HH-LH mixing 
occurs that can give rise to level anticrossings, 39 further in- 
creasing the complexity of holes' electronic properties. 

While in 2D systems the HH-LH mixing is absent at the 
subband edge and can be assumed to have a small magni- 
tude for finite k± in the presence of strong quantum con- 
finement, 50 this is not the case for hole nanowires 27 ' 5 1-55 
and quantum dots. 56-59 Previous studies of the topmost hole- 
nanowire eigenstates have shown that they are, to a varying 
degree, mixtures of HH and LH states. 55 ' 60 This has been 
found to influence, e.g., the spin splitting in cylindrical hole 
nanowires with hard-wall confinement 55 ' 60 ' 61 and nanosized 
clusters. 62 

In this article, we present an extensive theoretical study of 
the electronic structure of hole nanowires, focusing especially 
on the importance of HH-LH mixing and splitting. Further- 
more, we have investigated the ramifications of the spin-3/2 
nature of holes for the spin splitting of hole nanowire states, in 
particular elucidating the interplay between the bulk-material 
and bound-state orbital effects due to a magnetic field. 

We use the Luttinger description 47 of the semiconductor 
valence band to study the subband edges (k z = 0) of hole 
nanowires oriented along the z direction. In order to quanti- 
tatively analyze the strength of HH-LH mixing at the subband 
edge states, we make use of scalar invariants obtainable from 
a decomposition of the spin-3/2 density matrix in terms of 
a multipole expansion. 63 The squared dipole moment, which 
apart from a prefactor is a measure of the squared magnitude 
of the hole-spin polarization, is by its definition 63 independent 
of the quantization axis of the angular momentum used in the 
basis functions of the hole subband edge state. Thus, it pro- 
vides a rigorous, basis-independent way of characterizing the 
hole-spin polarization character of hole subband edge states. 

We present a comprehensive analysis of the hole-spin po- 
larization character of hole subband edge states in nanowires 
with different confining potentials. The strong HH-LH mixing 
manifests as large spatial variations in the hole-spin polariza- 
tion density profile of individual subband states. Furthermore, 
we observe large differences between the hole-spin polariza- 
tion density profiles of different subband states. We also find 
that the topmost hole subband edges in general display uni- 
versal features that are robust against changes of the confin- 
ing potential's symmetry. For hole subband edge states that 
are far away from the bulk valence band edge, however, dif- 
ferences in the strength and shape of the confining potential 
significantly affect the wavefunction's hole-spin polarization. 

A direct manifestation of the strong HH-LH mixing (and 
splitting) in hole nanowires is observed in the spin splitting 
due to a finite magnetic field. We use perturbative and nu- 



merical calculations to show that, in a magnetic field applied 
parallel to the wire axis, the effective g factor varies strongly 
between different wire subband edges. We analyze the sep- 
arate contributions to the g factor arising from bulk-material 
and orbital bound-state interactions with the magnetic field. 
We find that the hole envelope functions acquire large or- 
bital momentum due to the imposed quantum confinement, 
which strongly affects the effective g factor of hole nanowire 
states. Interestingly, the bulk-material and bound-state con- 
tributions can come with the same or opposite sign, thus ef- 
fectively enhancing or suppressing each other's contribution 
to the total spin splitting. We also find that the interplay be- 
tween the two contributions is very sensitive to the hole wave 
function, and thus, that the g factor of individual wire states 
can be tuned by changing the confining potential. For ex- 
ample, we show that by changing the aspect ratio of the lat- 
eral dimensions of a wire with rectangular crossection, one 
can tune the magnitude and even the sign of the g factor 
quite sensitively. Such a situation can, in principle, be eas- 
ily realized by tuning the lateral confinement via side gates in 
self-assembled 64 ' 65 or lithographically defined nanowires, 66 or 
quantum point contacts. 67-69 Quasi-onedimensional (ID) hole 
systems thus provide an interesting laboratory for the study of 
spin-3/2 physics, as well as have the potential for being build- 
ing blocks in hole-spintronics applications. 

In the following, we start by presenting the basic theoret- 
ical formalism for describing hole-nanowire states in the ab- 
sence of a magnetic field within the Luttinger model (Sec. II). 
Subsequently, we will discuss the influence of valence-band 
mixing on hole-nanowire subband edges (k z = 0) in Sec. Ill, 
using scalar invariants from the decomposition of the spin- 
3/2 density matrix to provide a comprehensive analysis of the 
hole-spin polarization of nanowire subband edges. A discus- 
sion of the influence of symmetry and strength of the confin- 
ing potential will be given. Next, the spin splitting of hole- 
nanowire subband edges due to an applied magnetic field is 
discussed in Sec. IV. We provide the theoretical descrip- 
tion of bulk-material and bound-state orbital interactions with 
a magnetic field in Sec. IV A, followed by a detailed anal- 
ysis of their individual contribution to the total g factor of 
hole-nanowire states in Sees. IV B and IV C. The effects of 
lower-symmetry corrections and finite k z will be considered 
in Sees. IV D and IV F. We discuss the important interplay 
between bulk-material and orbital contributions to the total g 
factors of hole nanowire states in Sec. IV E, demonstrating a 
versatile tunability of the magnitude and sign of the g factor 
by simple confinement engineering. We provide a summary 
and conclusions in Sec. V. 



II. THEORETICAL DESCRIPTION OF HOLE NANOWIRE 
SUBBANDS IN THE ABSENCE OF A MAGNETIC FIELD 

A. Luttinger model for bulk semiconductor valence bands 

The electronic structure of bulk valence bands can be de- 
scribed using the Luttinger Hamiltonian 47 
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where mo is the vacuum electron mass, fc, are the components 
of linear orbital momentum, Jj are the Cartesian angular- 
momentum operators for a particle with spin |, and ji the 
Luttinger parameters. In Eq. (1), we have used the notation 
{A, B} = \{AB + BA). We will choose our quantization 
axis for total (spin-3/2) angular momentum to be along the 
z-direction. In this representation, J z will be diagonal with 
eigenvalues ± | and ± \ . The majority of the results presented 
in this paper focus on the properties of subband edge proper- 
ties at k z — and, hence, we restrict our discussion of the 
theoretical formalism to this case in the following unless oth- 
erwise indicated. At the end of the paper, we will discuss 



hole-nanowire dispersions as a function of k z , as well as ad- 
dress the issue 52 ' 53 of subband edges occurring at finite k z . 

To provide clearer interpretations of our results, we will ini- 
tially use the Luttinger Hamiltonian in the spherical approx- 
imation 70 where band warping terms due to cubic-symmetry 
corrections 71 are neglected. Formally, the spherical approx- 
imation consists of setting 72 — > "f s and 73 — > 7 S with 
7 S = (272 + 373)/5. At the end of the paper, we will ex- 
plicitly discuss the effects of cubic corrections. In the spheri- 
cal approximation (and for k z = 0), the bulk Hamiltonian (1) 
simplifies to 
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where k± — k x ± ik y , k\ = k x + ky, and j± = 
(^J x ±iJ y ^J I \/2. The parameter 7 S is a measure of the 

strength of the SO coupling that is responsible for lifting the 
four-fold degeneracy for k ^ 0. This arises through the 
second term in Eq. (2), which describes the energy splitting 
between HH and LH bands. The third term is responsible 
for HH-LH mixing through the operators J±, which couple 
HH and LH amplitudes of the hole wave function. Thus, for 
k± 7^ and finite 7 S , the bulk valence bands can no longer be 
identified as being of purely HH or LH character in the sense 
of having definite spin projections J z = ±| or J z — ±\. 55 In 
the following discussions, we will restrict the use of the HH 
and LH labels solely for states that are pure HH or LH states 
in this sense. In the presence of a two-dimensional confin- 
ing potential, such as in a quantum wire along the z direction, 
even states at the band edges have a finite k± due to quantum 
confinement and thus are HH-LH mixtures. 51 ' 55 



lope functions is a constant of the motion. 70 71 This conserved 
operator F = J + L, which is similar to a total angular mo- 
mentum, was first used by Baldereschi and Lipari 70 ' 71 to sim- 
plify the acceptor-state problem utilizing an analogy between 
L and J and the L-S coupling scheme in atomic physics. Ser- 
cel and Vahala 72 ' 73 extended this approach to study spherical 
quantum dots and cylindrical wires. Their formalism will be 
the starting point for our own investigations of the properties 
of hole nanowires. 

We take the quantum-wire axis to be parallel to the spin- 
3/2 quantization axis (i.e., the z direction). The wire Hamilto- 
nian HI + V c h(r) [including the cylindrical hard-wall poten- 
tial: V c h(r) — for r < R and V c h(r) — 00 elsewhere] has 
common eigenstates with F z = J z + L z . Furthermore, it is 
a property of the Luttinger Hamiltonian that, at k z = 0, each 
block labeled by a fixed quantum number F z further decou- 
ples into two 2x2 blocks according to 



B. Description of hole states in nanowires 

We will now provide a theoretical description for sub- 
bands in hole nanowires confined by potentials with different 
strengths and confining potentials. We will start off by study- 
ing cylindrical hole nanowires confined by an infinite, hard- 
wall potential of radius R, which is particularly interesting 
due to its symmetry. We note that such wires can be readily 
fabricated by self-assembly. 18 

In the spherical approximation, the sum of the angular mo- 
menta J of the band-edge Bloch functions and L of the enve- 
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where the Hp are 2x2 matrices acting in the Hilbert sub- 
spaces spanned by states with spin projections J z = (§, — 5) 



(for a = +) and J z = (|, — |) (for a = -), respectively. 73 

To determine bulk-hole states at zero magnetic field (B z = 
0) that can be used to construct cylindrical-wire subband 
states, we use polar coordinates (r, <j>) and the wave function 
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ansatz 55 ' 73 



The corresponding eigenenergies are 
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where J„(fcj_r) is an integer Bessel function and 
are constants. Diagonalization of {H s L ) Fz in this representa 
tion yields the following bulk-hole eigenstates 
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In Eqs. (5) and (6) the subscript index v = ± labels eigen- 
states within the 2 x 2 subspaces spanned by HH and LH states 
with spin projections J z = ±| and J z — respectively. 
We emphasize that, for k± ^ 0, hole bands are no longer of a 
purely HH or LH type, hence we use the generic label v = ±. 
The coefficients , b^ , which do not depend on F z (nor 



on a at B, 



0), are given by 73 
1 



,(<*) = _. 
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h 2 k 2 
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A generalization of these results for finite z component of 
the magnetic field was presented in Ref. 55. There it was 
pointed out that an interesting crossover occurs with increas- 
ing k± for the effective g factor of bulk hole valence bands 
in the presence of a finite B z (see Fig. 2 of Ref. 55). For 
k± = 0, the spin splitting of the bulk-hole bands is character- 
ized (in absolute terms) by the g factors g = 6k and g = 2k, 
which follows from the form of the bulk-hole Zeeman Hamil- 
tonian that will be introduced in Eq. (15). However, at large 
k±, spin splitting turns out to be characterized by g = and 
g = 4k. This situation is reminiscent of the Zeeman effect 
in two-dimensional hole systems subject to an in-plane mag- 
netic field, 39 a result that will also be significant for our later 
discussion of the spin splitting of a particular class of cylin- 
drical hole-nanowire subband edges. 

Quantum wire eigenstates with subband index a are formed 
by superimposing bulk-hole eigenstates according to 



(9) 



The subband eigenenergies £^ ch ' and expansion coefficients 
c^ F v are obtained by imposing the hard-wall boundary con- 

(7) dition v P^ h ' (-R, 4>) = and solving the secular equations 



a 



(ch),(ch) 



br'J Fz - 3/2 (ktR)J Fz+1/2 (klR) 



2 (k-R)J Fz+1/2 (k+R) 

bf ) a^ ) J F ^ 2 {k+R)J Fz+m {klR)-b i ^ ) af\ 



>h),(ch) j 

a_ b + J Fz -3/ 

(<*) (ch) 



(10) 
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Again, a generalization for the case of finite magnetic field 
was presented earlier in Ref. 55. Hole nanowire subband 
edges can thus be obtained in a semi-analytical fashion for 
a cylindrical hard-wall confining potential. 

We have also investigated hole nanowires defined by con- 
fining potentials with other symmetry and strength, in order 
to gauge the universality of our findings. In particular, we 
have focused on two additional cases, hard-wall confined hole 
wires with square (and rectangular) crossection, and cylindri- 
cal hole nanowires confined by a harmonic potential. 

The first case corresponds to a confining potential that is 
given by 

V A (x, y) = for < x < W x and < y < W v 
V r h(x,y) = oo elsewhere , 



where W XtV are the lateral wire dimensions and W x — A ■ W y . 
We will first consider a square-crossection case with A = 1, 
deferring a discussion of the dependence on aspect ratio A 
to the end of our paper. Using a complete two-dimensional 
square-well basis set, the wave function of each individual 
wire subband edge with index a is expanded according to 



¥ a ^(x,y) = J2 



„(<*,&) 



J«,rh) 
\ j(<*,rh) 



w v 



(12) 



from which the hole wire eigenenergies _E( Q rh ) and expansion 
coefficients in Eq. (12) are obtained by numerical diagonal- 
ization of HI + V T h(x, y) within this basis. 
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Similar calculations are performed for cylindrical wires where uj is a measure for the softness of the harmonic poten- 
confined by a harmonic potential tial. The eigenstates of these wires are expanded in a complete 



basis set using Hermite polynomials according to 
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Numerical diagonalization of + V£ s (x, y) within this basis 
yields the corresponding eigenenergies 

£(«,cs) and 

expansion 

coefficients in Eq. (14). 

III. VALENCE BAND MIXING IN HOLE NANOWIRES 
A. Cylindrical hole wires with hard-wall confinement 

We will start by analyzing the hole subband edge states in 
cylindrical hard-wall wires, which are of the form given in 
Eq. (9). To begin with, we turn our attention to the detailed 
form of the secular equations (10) and (11). A special case 
arises for subbands with F z = o\. For these values of F z , 
the spinor entries in the hole-wire wave functions [Eq. (9)] 
are proportional to J±i(kj_r). Since J_i(x) = — Ji(x), the 
secular equations reduce to Ji(fc^r) = 0; a condition that can 
be satisfied by the individual bulk-hole eigenstates with v = 
± that are part of the superposition in Eq. (9). Thus, for F z = 
ah, the wire eigenstates are actually pure bulk states with v = 
+ or v = — . Note, however, that these bulk states are still 
of mixed HH and LH character for finite (quantized) fc^, 55 
hence our labeling in terms of generic quantum number v and 
not in terms of HH and LH as some authors do, reserving 
the latter labeling strictly for states that have definite angular- 
momentum projections J z = ±| and J z = ±~, respectively. 
For any other values F z = —<j\ and \F Z \ > ~ the indices 
of the Bessel functions appearing in Eqs. (10) and (11) are 
not equal and the resulting wire states are true mixtures of 
the bulk-hole eigenstates tp a F z + { r r, 4>) and t/j ff F ti -(r, <f>), for a 
given a, in order to satisfy the boundary condition set by the 
wire confinement. 

To clarify how strong the HH-LH mixing in fact is, we will 
analyze the wavefunction of hole nanowire subband edges in 
greater detail. To be specific, we use the Luttinger parameters 
characteristic for GaAs 74 in our calculations: 71 = 6.98,72 = 
2.06, 73 = 2.93, hence 7 S = 2.58. (We remind the reader 
that the following results are calculated using the spherical 
approximation for the Luttinger Hamiltonian.). At B z = 0, 
the wire eigenstates are two-fold degenerate, the degeneracy 
occuring between states corresponding to subspaces a = ± 
characterized by aF z . We will assume a = + in the following 



and suppress the index a for brevity. 

We illustrate the HH-LH mixing in two different ways. The 
traditional approach is to show the squared amplitudes of the 
HH and LH components in the subband eigenfunctions. How- 
ever, the thus obtained HH or LH density profile depends on 
the chosen basis, i.e., the quantization axis of total angular 
momentum. As an alternative, we will also consider a quan- 
tity that measures the invariant hole spin-polarization density, 
P = Pi/pq. Here the quantities p\ and p\ are scalar invariants 
derived from the spin-| density matrix 63 that are related to the 
hole spin-dipole density and total charge density, respectively. 
A pure HH state has P = 1.8 everywhere, whereas a pure LH 
character yields P = 0.2 uniformly in space. The formulation 
in terms of scalar invariants is basis-set independent. In prin- 
ciple, it therefore enables a more rigorous analysis of HH-LH 
mixing than would be possible in the traditional approach. 

In Fig. 2, results are shown for the topmost subband edge 
(i.e., the subband with index a = 1). The subband-edge 
state is characterized by F z = — |. Figures 2(a) and (b) give 
the cross-sectional profiles of the HH and LH components in 
the subband-edge wavefunction, illustrating that (i) the spatial 
density profile of the two components is very different, and 




FIG. 2: (Color online) Analysis of the highest-in-energy quantum 
wire subband-edge state with a — 1 and F z — — \ (i.e., a — +). 
Square amplitude of (a) the HH component and (b) the LH com- 
ponent of the subband edge wave function, (c) Total wavefunction 
amplitude squared, (d) Spatial variation of the hole spin-polarization 
density, P = pf/po- P = 0.2 corresponds to LH and P = 1.8 to 
HH character. 
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FIG. 3: (Color online) Analysis of the quantum wire subband edge 
with a = 2 and F z = i of a cylindrical, hard-wall confined hole 
nanowire. Squared amplitude of the (a) HH component, and (b) LH 
component of the subband edge wave function, (c) Total wavefunc- 
tion amplitude squared. The hole spin-polarization density, P, is zero 
for this subband throughout the entire crossection of the nanowire, 
and is thus not shown. 



(ii) their respective overall magnitudes differ by two orders of 
magnitude, with the LH component being the dominant one. 
These findings are also reflected in the (normalized) hole-spin 
polarization density, P, which is shown in Fig. 2(d). The nor- 
malization by /3q (which is a measure of the hole charge den- 
sity) ensures that the observed spatial variations in the hole- 
spin polarization density can indeed be attributed to hole spin 
and not charge-density variations. For comparison, the charge 
density profile is shown in Fig. 2(c). 

Inspection of the hole-spin polarization density P shows 
that the eigenstate is predominantly LH-like throughout most 
of the core of the wire, where P — 0.2. However, toward the 
edge, i.e., with increasing r, the hole-spin polarization van- 
ishes, only to recover to a finite value at the very edge of the 
cylindrical confining potential. This is a manifestation of HH- 
LH mixing. The predominantly LH character of the topmost 
hole-wire subband in GaAs is in agreement with other theo- 
retical predictions. 27 ' 54 ' 75 

We now turn our attention to the second-highest subband, 
with index a = 2. This wire subband edge is characterized by 
F z — | and therefore is actually a bulk-hole eigenstate. (See 
discussion above.) Figures 3(a)-3(c) show the squared am- 
plitudes of the HH and LH components of the wavefunction, 
along with the total hole charge density profile. The hole- 
spin polarization density turns out to vanish identically across 
the entire crossection of the wire and is therefore not shown. 
This means that the second-highest subband edge completely 
lacks spin polarization. Here the difference between the tradi- 
tional approach toward discussing HH-LH mixing and the al- 
ternative one based on spin-density-matrix invariants 63 is most 
transparent. We will elaborate on this fact more later on. 

The third-highest subband edge turns out to be of predom- 
inantly HH character, as can be seen from Fig. 4. The hole- 




FIG. 4: (Color online) Analysis of the quantum wire subband with 
a = 3 and F z = |. Squared amplitude of the (a) HH component, 
and (b) LH component of the total subband edge wave function, (c) 
Total wavefunction amplitude squared, (d) Spatial distribution of the 
hole spin-polarization density, P — pi/po- 



spin polarization density profile, Fig. 4(d), shows clearly that 
the wavefunction around the core of the wire is of purely HH 
character (P — 1.8). With increasing r, however, the hole- 
spin polarization drops to zero, subsequently recovering val- 
ues around P = 0.2, which signifies a LH character. The HH- 
LH mixture variations can also be seen explicitly by compar- 
ing the magnitudes and spatial profile of the square amplitudes 
of the HH and LH components of the wavefunction, shown in 
Figs. 4(a) and (b), respectively. 

For subbands at decreasing energy (i.e., those with higher 
subband index a), the hole-spin polarization density profiles 
become increasingly complex. In Fig. 5 we show the hole-spin 
polarization density P for the subband edges with a = 4 ... 8 
and a = 10. The subband edge with a = 9 has F z = | 
similarly to the one with a = 2 and, hence, is again also 
a bulk-hole eigenstate with vanishing polarization across the 
entire wire crossection. 

The fourth-highest level (a = 4) is predominantly LH-like 




x/R x/R x/R 



FIG. 5: (Color online) Spatial distribution of the hole spin- 
polarization density, P = pi/pi, of subband edge states with (a) 
a = 4, F z = -§ (b) a = 5, F z = -§ (c) a = 6, F z = § (d) 
a = 7, F z = | (e) a = 8, F z = -f (f) a = 10, F z = -§. The 
wire subband with a — 9 is characterized by F z — 5 and has zero 
hole-spin polarization throughout the entire crossection of the wire, 
and is thus not shown. 
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FIG. 6: (Color online) Spatial distribution of the hole spin-polarization density, P — pi/po for the levels a = 1,2,3,6,8, for square- 
crossection hard-wall confinement (a)-(e) and cylindrical soft-wall confinement (f)-(j) hole nanowires. R a — [2E^ a ' cs ^ / (■yi + 5/2j a )] 1 ^ 2 . 



throughout the core of the wire. The hole-spin polarization 
drops to zero with increasing r, only to recover and switch to 
a HH character, P = 1.8, at the edge of the wire. A similar 
behavior is observed for the subband edge with a = 8. The 
other subband edges show equally large variations of the hole- 
spin polarization as a function of the radial coordinate r. This 
has been found to be a general feature for all subband-edge 
levels characterized by F z ^ ~ and is a clear manifestation of 
the strong valence band mixing. 



B. Dependence on shape and strength of the wire confinement 

To investigate the universality of our findings, we have per- 
formed similar analysis for hole nanowires with confining po- 
tentials of different shape (or symmetry) and strength. In par- 
ticular, we considered two additional cases, namely hard-wall 
confined hole wires with square crossection, and cylindrical 
hole nanowires confined by a harmonic potential. We con- 
tinue using the spherical approximation for the moment, de- 
ferring a discussion about cubic crystal-symmetry effects to 
later parts of the paper. 

The subband edges of wires defined by the two additional 
types of confinement are analyzed in the same way as we 
did previously for the cylindrical, hard-wall confined ones. 
In particular, we are again focusing on the hole-spin polar- 
ization density profiles to reveal the nature of HH-LH mix- 
ing. In Fig. 6, we show the calculated hole-spin polariza- 
tion density profiles, P, for square crossection, hard-wall (up- 
per row), and cylindrical, soft-wall hole nanowires (bottom 
row). For brevity, we only focus on the subband edges with 
a = 1, 2, 3, 6, 8 to illustrate our main findings. 

The topmost subband with a = 1 is seen to be of LH 
character, i.e., having P = 0.2, throughout most of the wire 
crossection. Especially the cylindrical wire with harmonic 
confining potential shows a hole-spin polarization density pro- 
file that is very similar to the corresponding hard-wall, cylin- 
drical wire. This is not surprising, as both confining potentials 
are axially symmetric. For the square-crossection wire, on the 



other hand, the lower symmetry of the confining potential is 
reflected in the polarization density profile. Nevertheless, the 
predominant character of the subband edge is still LH-like, 
showing that this feature is indeed very robust. 27 ' 54 

The second-highest subband edge (with a = 2) shows 
quite a different hole-spin polarization profile when compar- 
ing the different confining potentials. The square crossec- 
tion nanowire shows a vanishing polarization at the center 
of the wire, similar to the cylindrical hard-wall wire. How- 
ever, due the reduced symmetry, the polarization is finite at 
the corners of the wire, where the hole-spin polarization is 
of LH-character. For the soft-wall cylindrical wire, the sec- 
ond subband displays a predominantly HH character through- 
out most of the nanowire crossection. However, examining 
the third-highest subband edges with a — 3, we see that 
the square crossection, hard-wall wire state exhibits a HH- 
character (similar to the cylindrical hard-wall case), whereas 
the cylindrical soft- wall nanowire has zero polarization. Thus, 
a level switching in energy between the subband edges with 
a = 2 and a — 3 has occured for the soft-wall cylindri- 
cal nanowire in comparison with its hard-wall cylindrical and 
square crossection counterparts. 

For subband edges with higher a, the deviations from ax- 
ial symmetry for the square crossection wire become increas- 
ingly more prominent, as is shown explicitly for the subband 
edges with a = 6 and a — 8. The soft-wall cylindrical wire 
on the other hand displays hole-spin polarization profiles sim- 
ilar to the hard-wall cylindrical case. 

Quite generally, our extensive studies of the hole-spin po- 
larization of subband edges in nanowires show that HH-LH 
mixing strongly affects the polarization of individual sub- 
bands, a characteristic that is universal and independent of 
the symmetry and shape of the confining potential. In addi- 
tion, we find that the topmost subband edge states with small 
a have HH-LH characters that are robust against changes in 
the confining potentials. 
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IV. SPIN SPLITTING OF HOLE-WIRE SUBBAND EDGES 
A. Luttinger model for bulk holes subject to a magnetic field 

We will now focus on the spin splitting of hole nanowire 
states in the presence of a magnetic field applied parallel to 
the wire axis. The external field couples to both the spin-3/2 
and orbital degrees of freedom. The Zeeman Hamiltonian de- 
scribing the interaction between the magnetic field, B z , and 
the spin degree of freedom is given by 

H z = -2k^i b B z ■ J z , (15) 

where ps is the Bohr magneton and n is the bulk hole g fac- 
tor. 39 We will neglect the anisotropic Zeeman contribution 



that is of higher order in J because, in typical semiconduc- 
tors and for the crystallographic wire direction we consider, it 
is much smaller than the one given in Eq. (15). 

In order to include the orbital effects due to a magnetic field 
applied in the z direction, we replace k — > k + eA. The 
symmetric gauge A = (— f , §,0)B Z will be used, and we 
will only consider resulting terms linear in B z because we 
extract g factors from the small-magnetic-field limit. In the 
spherical approximation, i.e., using the Luttinger Hamiltonian 
of Eq. (2), the orbital terms arising from the replacement of 
canonical wavevector with the kinetic one have the following 
form (we use atomic units and the definitions L z = xk y — yk x 
and x± — x ± iy) 



H, 



orb 



H 



orb.diag 



H, 



orb, mix 



71 • l 4 x4 +7s \J Z - -■ 14x4 



L z ^bB z - i"f s (x-k-Jl - x+k+J 2 ^ fJ,sB z 



(16) 



The first term is proportional to the orbital angular momentum 
operator L z and is diagonal in spin space. We denote this con- 
tribution by H^ lh diag . The remaining terms are off-diagonal 
in spin space, coupling states with J z eigenvalue differing by 
2. We denote this contribution by H^ lh mix . These two contri- 
butions will be analyzed in greater detail in our discussion of 
the orbital contribution to the total g factor of hole nanowire 
subband edges. 

At zero magnetic field, subbands with a given a and uF z 
are doubly degenerate. A finite B z lifts this degeneracy, giv- 
ing rise to a spin splitting from which we extract a g factor 
according to 



H. 



orb, i 



and H* vh 



9to 



B lirn o Sgn((F,) + ) 



E+{B Z ) - E~{B Z ) 
PbB z 



(17) 



The sign factor introduces the physical definition of assign- 
ing sign according to the projection of (F) along the magnetic 
field direction. It removes any ambiguity relating to the par- 
ticular choice of symmetry, a, labelling the 2x2 blocks of 
the Hamiltonian. Due to the finite energy separation between 
wire subbands, this definition is equivalent to the perturbative 
(in B z ) result 

n a - gq-nrfV F \ + ^ ^ Z + g orb)« ~ ( g Z + H^ h )~ 
i/tot ^ °B n \\ r z/ a ) HbB 1 



ts . \ + 



sgn((F 2 ) + ) 2(gz+ ^ b/ " 

PB&z 



(18) 



where the expectation values are taken with respect to the 
B z = wire subband edge eigenstates with a = +. 

In the following, we analyze the individual bulk-material 
and orbital contributions to the total g factor gf ot — g% + 
2orb,dia g + 5orb,mix that arise from the linear-in-S z terms H z , 



, respectively. It is useful to under- 
stand each contribution's separate properties, as their relative 
importance may vary depending on physical details. For ex- 
ample, the bulk-material contribution will dominate in di- 
luted magnetic semiconductor wires where k can be enhanced 
by several orders of magnitude due to the p-d exchange inter- 
action between hole carriers and magnetic dopant ions. 60 ' 76 



B. Bulk-material contribution to g factors 

We begin by analyzing the contribution to g™ ot that is due to 
the Hamiltonian Hz- It follows from Eqs. (18) and (15) that 
the g factor derived from the bulk-hole Zeeman effect has the 
form 

g z - -4«sgn((F z } + ) = -2V$k sgn((F z )+)p 10 . 

(19) 

In the r.h.s equality, we have expressed g% in terms of pio, 
which is a dipole-related component of the spin-3/2 density 
matrix (defined in Table III of Ref. 63). 

It is clear from Eq. (19) that an eigenstate of pure HH 
character has an effective g factor of —6k, while a pure LH 
state has a spin splitting characterized by g = —2k. This 
is certainly the case for bulk-hole states at the valence-band 
edges. 39 The situation is different for quantum-confined hole 
states, e.g., in a wire. 5155 As discussed in Sec. Ill, strong HH- 
LH mixing is a universal characteristic of hole nanowire sub- 
band edges, and hole nanowire states of purely HH or LH 
character are not expected to be found. In particular, the anal- 
ysis of the squared amplitude of the spin-3/2 dipole moment 
in Sec. Ill showed that the normalized hole-spin polarization 
profiles, P = Pi/po display strong spatial variations as a man- 
ifestation of HH-LH mixing, (pi is directly related to p w ; see 
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FIG. 7: (Color online) Bulk-material contribution to the (/-factor of 
the highest-in-energy subband-edge states, with a = 1 ... 10, of a 
hole nanowire defined by a cylindrical hard-wall confinement. The 
arrows indicate the values of gz/ft = —2 and —6 corresponding to 
the bulk-hole (/-factor of pure LH and HH states, respectively. 



Ref. 63 for details.) The question is thus: To what extent does 
the HH-LH mixing change the g factors for hole states in a 
wire as compared with the bulk-material (HH and LH) values 
of —6k and —2k? 

Figure 7 shows the effective g factors normalized by k, of 
the ten highest subband edges in cylindrical nanowires with 
hard-wall confinement. The values display strong variations 
as a function of subband index a. Furthermore, most sub- 
bands have values that are strongly suppressed in comparison 
to the values expected for pure HH and LH states, respec- 
tively. Inspection of the hole-spin density profiles discussed 
in Sec. Ill and comparison with the (/-factor values shows a 
strong correlation, establishing a direct connection between 
HH-LH mixing and the observed large (/-factor variations as a 
function of subband index. For example, the topmost subband 
edge (a = 1) was found to be predominantly of LH character 
(see Fig. 2). Correspondingly, its g factor is found to be close 
to —2k, just as expected from Eq. (19) for a pure LH state. 

The second (and ninth) subband edges from the top have 
vanishing polarizations, which correlates with a vanishing g 
factor. Both the second and the ninth state belong to the class 
of subband edges with F z — \ and are thus also a bulk-hole 
eigenstate. As discussed in Fig. 2 of Ref. 55, the spin splitting 
of bulk hole bands is characterized at high values of k± by 
values of g = or \g\ = 4k. The subbands with a = 2 and 
a = 9 are examples of the former type. We have also found 
wire-subband edges of the other type whose spin splitting is 
characterized by \g\ = 4k, e.g., the one with a = 16. 

The dipole moment turns out to be just one of three non- 
trivial multipole moments that characterize a spin-3/2 state. 63 
This is a fundamental difference with a spin- 1/2 system where 
the dipole moment suffices to uniquely determine the spin 
state. When the spin-3/2 polarization (i.e., the dipole moment) 
vanishes for pure spin states in two-dimensional 63 and one- 
dimensional 55 hole systems, a substantial octupole moment 
can exist. For hole nanowire subband edges, a mathematical 



TABLE I: Bulk-material contribution to g-factors for hole nanowire 
subband-edges with index a = 1 ... 10 for three different types of 
wire confinement: (i) hard-wall cylindrical cross-section, (ii) hard- 
wall square cross-section, and (iii) soft-wall cylindrical cross- 



section. 



a 


1 


2 


3 


4 


5 


6 


7 


8 9 


10 


0) 9z 


-2.26 


0.00 


-2.66 


-1.70 


0.19 


-4.13 


-0.65 


-1.20 0.00 


-0.28 


0i)sf 


-2.22 


0.07 


-3.39 


-1.55 


0.06 


-2.58 


-0.72 


-0.59 0.08 


-0.17 


(iii)ffS 


-1.84 


-5.57 


0.00 


-1.34 


-0.66 


-1.60 


-2.68 


-0.97 0.01 


-0.94 



relation illustrates this point: 



Pa = 2 P2 = Pi + Pi 



(20) 



Equation (20) links the scalar invariants pf obtained from 
the multipole-like expansion of the spin-3/2 density matrix. 63 
The squared monopole, p\, and dipole, p\, correspond (apart 
from prefactors) to the squared magnitude of the hole charge 
density and hole-spin polarization density, respectively. The 
quadrupole and octupole are unique to spin-3/2 systems. The 
quadrupole is a measure of the HH-LH mixing, while the oc- 
tupole does not have a straightforward physical interpretation. 
The left equality in Eq. (20) quantifies the HH-LH mixing that 
is present at subband edges. The right-hand equality implies 
that the hole-spin polarization, measured by the dipole mo- 
ment pi, and the octupole, p$, are complementary. Thus, for 
the subband edges for which the hole-spin polarization van- 
ishes, a maximum octupole moment is induced. Similar com- 
parisons can be made for all the wire eigenstates. In general, 
states with mainly HH (LH) character will have g factors close 
to —6k (g — —2k). States with mixed character, displaying 
fluctuations in the hole-spin polarization profiles, have widely 
varying values for g. 

In Table I the calculated values of the g factors (with 
k = 1.2, the value 39 corresponding to GaAs) for hole wires 
with square-crossection hard- wall confinement (mid-row) and 
cylindrical harmonic potential (bottom row) are shown along 
with the numerical values for the (/-factors corresponding to 
Fig. 7 for the cylindrical hard-wall nanowire (top row). In 
general, we observe strong similarities between the g-factors 
of the different structures, with respect to both magnitude and 
sign. This is especially clear at low energies, where the axial 
symmetry seems to be retained even in the wire with square 
crossection. These findings do not come as a surprise when 
one considers the similarity of hole-spin polarization density 
profiles for the low-lying subband edges (shown in Fig. 6). As 
could be expected, the similarity is stronger for the two types 
of cylindrical confinement, while the square-crossection wire 
states display some complex hole-spin polarization variations 
(e.g., for subband edges with a — 6 and 8), with correspond- 
ing discrepancies between the calculated g factor values of 
the square crossection wire and the corresponding cylindrical 
ones. Note also that the subband edges with vanishing polar- 
ization in the cylindrical wires, acquire a finite, albeit small, 
polarization for the square crossection one, with a correspond- 
ing non-zero value of the hole g factor. Another interesting 
observation is the level switching that occurs between sub- 
band edges with a — 2 and 3 for the cylindrical wire with 
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TABLE II: Total orbital g factors, g° b , and the individual contributions gSb,dia g ar >d gSb,mix f° r tne ten highest subband edges with a = 1 ... 10 
in cylindrical hardwall hole nanowires. 

q,gFjl,-§ 2,\ 3,f 4,-f 5,-| 6,| 7,| 8,-§ 9,f 10,— § 

2.78 -1.82 4.91 -0.05 -8.62 -14.10 -4.46 -3.28 -1.82 -18.09 

ffSbdiag -0-54 -1.82 -8.27 -12.40 -10.33 -5.66 -24.10 -25.04 -1.82 -19.63 

gZ bmix 3.32 0.00 13.18 12.35 1.71 -8.44 19.64 21.76 0.00 1.54 



parabolic confining potential. This is an example for the sen- 
sitivity of physical quantities (here the wire bound-state en- 
ergy) with respect to quantum-confinement-induced HH-LH 
mixing. 

We conclude that, as a universal characteristic of hole 
nanowires, states are of mixed HH and LH character even 
at the subband edges. This is manifested by large variations 
in the bulk-material contribution to g factors as a function of 
subband index, concomitant with an overall suppressed mag- 
nitude of Zeeman splitting as compared with free holes. 

C. Orbital effects 

We now turn our attention to how the coupling between the 
holes' orbital degrees of freedom and the magnetic field af- 
fects the splitting of hole-nanowire subband edges in a mag- 
netic field applied parallel to the wire axis. Using a similar 
perturbation approach as in Sec. IV B , the orbital contribu- 
tion to the total g-factor for a subband edge with index a can 
be written as 



where 



Sorb, diag 
i/orb,mix 



where the angular brackets (|||) and (— 1| — i) are squared 
amplitudes of the HH and LH components of an eigenstate 

ip+F^ as described in Eq. (9). Equation (25) is linear in F z , 
reflecting the angular momentum that is acquired by the enve- 
lope wave function of the hole subband edge due to the cylin- 
drically symmetric quantum confinement. We also note that 
the two terms in Eq. (25) differ by a factor of 71 ± j s . This is 
a manifestation of HH-LH splitting. 



The total orbital g-factor contribution g" rb is shown in Ta- 
ble II (top row) for the ten highest-in-energy subband edges 
(a = 1 . . . 10), of a cylindrical hard-wall wire. The orbital g 
factor displays large variations as a function of subband index 
a. Moreover, for some subband indices, e.g., a = 6 and 10, 
the values are significantly larger than the corresponding bulk- 
material contributions discussed in the previous subsection. In 
general, we observe that the orbital contribution to the total g 
factor of hole subband edges in cylindrical nanowires domi- 
nates over the corresponding bulk-material contributions. 

Examining the individual contributions go rbidiag and 
forb mix' we see tnat ' ^ rst °f tne two contributions have 
similar magnitudes. In addition, the two contributions have in 
general opposite sign, with the exception of that for the sub- 
band edge with a = 6. The overall sign of the orbital g" lh is 
thus determined by the competition between these two terms. 

The origin of large fluctuations in the total orbital contribu- 
tion to the g factor can be understood by examining the inter- 
play between the diagonal and mixing terms in greater detail 
for cylindrical hard- wall hole wires. The operators entering 
Eqs. (22) and (23) can be expressed in polar coordinates as 



I 

In Fig. 8(a), we show values of g" rb diag for subband edges 
with different F z . For each value of F z , the four g-factor val- 
ues correspond to the four highest-in-energy subband edges 
with that particular F z . We label this sequence in energy by 
the index (3, not to be confused with the general subband in- 
dex a. A linear dependence in F z can be clearly seen from 
the figure, confirming the presence of the large quantum- 
confinement-induced envelope angular momentum. We also 



-2 



Sorb — Sorb, diag + Sorb, mix ' 

L z i ix4 ) + a + Is ((4 2 - fl4x 4 ) L 



71 



sgn{(F z )+) 



sgn«F z >+) 
(«x_fc_J+ — ix + k + J 2 



&2) 



L z = -idf , 
£_A;_ = -iL 2 _{rdr + L z ) 
x + k + = (a;_fc_)^ , 



(24) 



Using the wire eigenstates shown in Eq. (9) it follows that, in 
the limit B z — * 0, the contribution of H^ rb diag to the orbital 
(23jf)factor is given by 



.9orb,dia g - -2sgn((F^) + ) 



F z -- 2 j (71+7 .)/-|- 



+ [F z + l) ( 7l - 7s )(-l|-l) + 



(25) 
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FIG. 8: (Color online) Values for individual orbital g-factor contribu- 
tions (a) fforb.diag and 0>) s£b,mix for subband edges with -| < F z < 
(for the cr = + states) in a cylindrical hard-wall hole nanowire. 
The index /3 labels subband edges with a given f z in order of de- 
creasing energy. 



observe a spread in the g-factor values for a fixed value of 
F z corresponding to different energies (and subband index a). 
This is a clear manifestation of HH-LH mixing, which enters 
through the squared amplitudes of the HH and LH compo- 
nents in Eq. (25). Thus, even though the Hamiltonian giving 
rise to the splitting characterized by <?" rb diag is diagonal in 
spin-space, its contribution to the spin splitting is affected by 
the quantum-confinement-induced HH-LH mixing in the un- 
perturbed states. The fluctuations in g" lh diag as a function of 
subband index a are therefore due to the different subbands 
having different total angular momentum F z , as well as dif- 
ferent HH-LH character. 

In contrast to H^. h diag , the Hamiltonian H^ lh mix couples 
hole subband edge wave function components with different 
J z . Thus, it gives rise to HH-LH mixing that is in addition 
to the HH-LH mixing inherent to the unperturbed wire eigen- 
states discussed in Sec. III. For a cylindrical hard-wall hole 
wire, this generates a contribution to the spin splitting given 

by 



5orb 



-4V3 



sgn«F 2 ) + ) \2 



(rd r 



L- 



(26) 

Here the notation indicates that the expectation value is be- 



tween HH and LH components of the subband edge state with 
a and +F Z , with spin projections J z = 3/2 and J z = —1/2, 
respectively. 

In Fig. 8(b), we show the g-factor contributions g" th mix 
for subband edges with — | < F z < Again, for each 
value of F z we show the calculated g factor values for the four 
highest-in-energy subband edges with the given F z . While the 
characteristics of the mixing term are more complicated, we 
can discern some general trends. For instance, for the lowest- 
in-energy subband edges, with (3 = 1, a clear linear depen- 
dence on F z is seen in the g factors. However, these g factors 
increase monotonically with \F Z \, in contrast to the splitting 
term, ff^bdiag' which 1S monotonically decreasing with \F Z \. 
In general, the terms fl£ rb)diag and g°. h ^ ix have similar mag- 
nitudes but opposite sign. However, for higher energies, this 
general rule does not apply and a large spread between the g 
factors for a given F z at different energies can be observed. 

A comparison between the orbital contribution to g factors 
fforb calculated for our three model hole nanowire structures 
is shown in Table III. Interestingly, with the exceptions of 
the levels with a = 2 and 3, all nanowires have very similar 
orbital g factors for low-index levels (small a). This is true 
both for the magnitude and the overall sign of g£ lh . In par- 
ticular, the square-crossection wire displays similar values to 
the cylindrical ones, despite the lower symmetry of the con- 
fining potential. Closer inspection of the levels with a = 2 
and 3 also shows that, in fact, the similarities extend to these 
levels as well, if one interchanges the values for a — 2 and 3 
for the cylindrical wire with harmonic confinement. This is in 
agreement with our previous discussion on the level switching 
observed in the hole-spin polarization of these levels for this 
wire confinement potential; see Sec. III. The spread between 
g" lh factors with given a corresponding to different confining 
potentials increases, however, with increasing subband index 
a, in accordance with the expectation that the more strongly 
delocalized wave functions probe the outer edges of the con- 
fining potential. 



D. Cubic corrections 

Next we consider the influence of corrections due to cubic 
crystal symmetry. These corrections depend on the orienta- 
tion of the wire with respect to crystallographic axes. In the 
following, we assume that the wire and, hence, the spin-3/2 
quantization axis, are parallel to the [001] direction. Starting 
from the full Luttinger Hamiltonian Hl shown in Eq. (1), we 
perform the transformation k — > k + eA in the symmetric 



TABLE III: Total orbital g factors, g" lb , for the ten highest subband 
edges (a — 1 . . . 10) in (i) cylindrical hard-wall, (ii) square hard- 
wall, and (iii) cylindrical soft-wall hole nanowires. 



a 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


0) 3orb 


2.78 


-1.82 


4.91 


-0.05 


-8.62 


-14.10 


-4.46 


-3.28 


-1.82 


-18.09 


(ii) Sort, 


2.74 


-2.21 


5.88 


-0.69 


-6.26 


-12.12 


-2.39 


-1.66 


-1.91 


-9.26 


(iii) So* 


3.39 


5.88 


-1.82 


-1.13 


-3.38 


-9.55 


-1.81 


-4.67 


-1.66 


-6.85 
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gauge. This yields the following terms linear in B z : 



H rb — — 



71 ' 14x4 + 72 I J; 



■ 72 



(ykx 



Ji+Jt)+i 



«73 (xk x - yky^j ( 



Jl - Ji 



r 



(27) 



In contrast to Eq. (16), the orbital Hamiltonian of Eq. (27) 
contains terms proportional to 72 and 73. The fact that 72 7^ 
73 accounts for bandwarping effects. 

In Fig. 9, we show the calculated values for g" ot for the 
hole wire subband edges with a — 1 ... 10 of a square cross- 
section hard-wall nanowire. The open symbols are values 
obtained from the numerical diagonalization of i/£ + Hz + 
_fforb + V rri (x,y), i.e., using the spherical approximation for 
the Luttinger Hamiltonian. The filled symbols are obtained by 
diagonalizing Hi, + Hz + H OID + V iri (x, y), which includes 
band warping due to the cubic crystal symmetry. For subband 
edges with low values of a, the g factors obtained from the 
two approaches turn out to be very similar. For higher val- 
ues of a, on the other hand, we observe an increasing spread 
between the g-factor values corresponding to the two cases. 
Properties of wires aligned parallel to crystallographic direc- 
tions with lower symmetry than [001] can be expected to be 
more strongly affected by band-warping corrections. 40 



E. Interplay of bulk-material and orbital wire-bound-state 
contributions to g factors: Tunability of sign and magnitude 

We now discuss in greater detail the interplay between 
the various contributions to the total nanowire-hole g factor. 



Lets recall that g" ot is the sum of three terms, g^ Q 



9z 



a , a 

?orb,diag ' yorb,mix 
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FIG. 9: Hole g factors g" ot for the top ten subband edges (a — 
1 . . . 10) of a square-crossection hard-wall nanowire. Circles cor- 
respond to values obtained from diagonalizing HI + H^ lh + Hz + 
V r h{x,y), i.e., using the spherical approximation for the Luttinger 
Hamiltonian. 70 Squares denote the values obtained from numerically 
diagonalizing Hl + H OI h + Hz + V t h(x, y), i.e., using the Luttinger 
Hamiltonian reflecting the cubic crystal symmetry. 



spin-space orbital, and HH-LH mixed orbital contributions, 
respectively. In Table IV, we show the calculated g" ot and 
its components for the ten highest subband edges in hole 
nanowires with square-crossection hard-wall quantum con- 
finement. The values are obtained by numerical diagonaliza- 
tion of Hi, + Horb + Hz + V T h(x, y), i.e., using the Luttinger 
Hamiltonian with cubic (not spherical) symmetry. 

Table IV illustrates several of our main conclusions. The 
bulk-material contribution shows strong variations in sign and 
magnitude as a function of subband index. This is a manifesta- 
tion of HH-LH mixing at hole-nanowire subband edges. Sim- 
ilar fluctuations occur for the orbital contributions. This is in 
part due to the fact that the envelope function of subbands with 
different a acquire different quantum-confinement-induced 
orbital bound-state angular momenta. In addition, the fluc- 
tuations also originate from HH-LH mixing. The two orbital 
contributions g" rb diag and fl& bjjnix have, in general, compa- 
rable magnitudes but opposite sign. 

Most importantly, the total g factor g" ot depends sensitively 
on the relative sign and magnitude of the three contributions. 
As both the orbital contribution to the spin splitting and the 
detailed form of the HH-LH mixing depend on the details 
of the quantum confinement, 56 we expect that the interplay 
between the bulk-material and orbital ^-factor terms can be 
tuned by confinement engineering. This would enable the use 
of nanowires to realize devices for nanospintronics. In Fig. 10 
we demonstrate such tunability of hole g factors in nanowires. 
The thick curves show g" ot for the three highest-in-energy 
subband edges of rectangular-crossection hard-wall confined 
wires plotted as a function of the aspect ratio A = W x /W y . 
The bulk-material and orbital contributions are also shown (by 
the thin solid and dashed lines, respectively). A general char- 
acteristic that is observed for all three subband edges is that 
the total g factor approaches zero in the limit of large aspect 
ratio A. This is due to the fact that, with increasing A, our 
nanowire system approaches the 2D limit, essentially mimick- 
ing a quantum well subject to an in-plane magnetic field. For 
such a system, the topmost hole-quantum-well subband edge 
is known to have a small (within the spherical approximation: 
vanishing) g factor, 39 ' 40 in agreement with our findings. The 
sudden sign change observed in g 2 and g 3 occurs at anticross- 
ing points in the energy spectrum where (F z ) goes through 
zero, and consequently reverses sign. 

The observed dependence of the subband edge g factor val- 
ues on the aspect ratio of the confining potential is quite in- 
triguing. All subband edges display a strong tunability of 
the magnitude of the g factor with a change in the aspect ra- 
tio. For the third subband, a variation of A between 1 and 2 
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TABLE IV: Total g factor g" ot — g% + gorb.diag + Sorb, mix f° r me ten highest-in-energy subband edges (q = 1 . . . 10) of a square-crossection 
hole nanowire with hard-wall confinement, is the bulk-material contribution, which arises from the Zeeman effect of holes in the semi- 
conductor material. The additional spin splitting due to coupling of the applied magnetic field to the orbital wire bound state is embodied in 
Sorb.diag an d Sorb, mix' which are derived from the first and second terms in Eq. (27), respectively. Values shown are obtained by numerically 
diagonalizing the cubic Hamiltonian Hl + // or b + Hz + V t h(x, y) and considering the spin splitting for B z — > 0. 
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2.61 
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17.52 


19.58 
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22.47 



gives rise to a change in the sign of the g factor that is not 
associated with an anticrossing, thus demonstrating the abil- 
ity to gradually tune the spin splitting to zero. The possibil- 
ity to manipulate both the magnitude and sign of g factors 
is a useful ingredient for spintronics applications. We there- 
fore propose that hole nanowires may be versatile building 




a Z 3 orb 
4 5 6 7 

A=W X /V\^ 

FIG. 10: Dependence of subband-edge g factors g" ot on the aspect 
ratio A for a rectangular-crossection hard-wall-confined nanowire 
(thick solid curves). We show results for the top three hole- wire sub- 
band edges: (a) a = 1, (b) a = 2, and (c) a — 3. The thin solid 
and dashed curves are plots of the corresponding bulk-material and 
orbital contributions, g^ and <7o lb , respectively. 



blocks for nanospintronics due to the demonstrated tunabil- 
ity that arises from the interplay between bulk-material and 
orbital spin splittings in an applied magnetic field. Confine- 
ment engineering of nanowires and quantum point contacts 
can be readily achieved with present-day technologies, e.g., 
using side gates on lithographically defined quantum wires 66 
and quantum point contacts, 41 ' 42 or wrap-gated self-assembled 



64.65 



nanowires. 



F. Hole nanowire dispersions for finite k z 

We finish our analysis by discussing the energy dispersion 
of hole nanowire states due to motion along the wire direction. 
For finite k z , the decomposition shown in Eq. (3) is no longer 
possible. Instead, the full 4x4 Luttinger Hamiltonian must 
be considered. We diagonalized numerically the Hamiltonian 
consisting of the sum of Hi shown in Eq. (1) and V v h(x, y) for 
a square-crossection hard-wall GaAs nanowire oriented along 
the [001] direction. 

The highest hole-nanowire subband dispersions are shown 
in Fig. 11. Several observations can be made. First, the 
hole dispersions are strongly non-parabolic, which is a distinct 
manifestation of the SO coupling and HH-LH mixing in the 
valence band of typical semiconductors. Second, some sub- 
bands display electron-like dispersions as well as off-center 
(k z 7^ 0) maxima. For example, the highlighted points A, 
B, and D correspond to off-center maxima for the subbands 
with a = 3,4 and 8, respectively. The occurrence of sub- 
band edges away from the zone center, as well as the strong 
non-parabolic dispersions, will have direct implications for 
the transport and optical properties of hole nanowires. 

One example is the quantized conductance of one- 
dimensional systems. It is well known that transport mea- 
surements in the linear-response regime of one-dimensional 
quantum wires and quantum point contacts yield quantized 
conductance steps 66 ' 77 78 if the transport is ballistic. Such con- 
ductance steps occur whenever a one-dimensional subband is 
opened as a channel for conduction, as it falls below the elec- 
trochemical potential determined by the contacts. Typically, 
top or side gate are employed to adjust the wire width and, 
thus, the quasi-lD subband energies. In the absence of a mag- 
netic field, the conductance steps occur in units of 2e 2 /h. 

In Fig. 11(b), we show a schematics of the linear-response 
conductance of a hole nanowire with the dispersion shown in 
Fig. 11(a). The abscissa refers to a gate voltage applied to a 
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FIG. 11: (a) Zero-field hole-nanowire subband dispersions for a 
square crossection hard-wall nanowire. The energy has been nor- 
malized by Ew = jih 2 /(2moW^). The hole nanowire states were 
obtained by numerical diagonalization of the cubic-symmetry Hamil- 
tonian Hi. + V r h{x,y). Dashed arrows indicate extremal values of 
subband dispersions that would be associated with steps in the two- 
terminal hole-wire conductance, (b) Schematic dependence on side- 
gate voltage for the linear two-terminal conductance through a hole 
nanowire or quantum point contact, based on the hole-subband dis- 
persions shown in (a). Each conductance step corresponds to the 
Fermi energy reaching a value indicated by the dashed arrows. 
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FIG. 12: (Color online) Absolute values of hole-wire g factors g" ot 
for states at finite wave vector k z for motion parallel to the wire axis. 
Results are shown for the subbands with a = 1 ... 5. 

between hole-like and electron-like dispersion, as is the case 
for point F. A measurement of the transport properties in the 
linear-response regime can thus yield insight to the detailed 
structure of the hole nanowire subband dispersions. 

In the presence of a magnetic field, the degeneracy of the 
hole nanowire subbands is lifted. Similar to the so far dis- 
cussed spin splitting of the subband edges at k z — 0, the spin 
splitting of the off-center maxima can also be probed by trans- 
port or optical measurements. For the square-crossection wire 
with dispersion shown in Fig. 1 1(a), we have three off-center 
maxima; at points A, B, and D. Calculation of the spin split- 
ting at these three points for a small magnetic field yields the g 
factor magnitudes gf oi = 2.25, gg t = 2.89, and g° t = 4.20, 
respectively. The full dependence of the g factors of states 
in the five highest hole-wire subbands as a function of k z is 
shown in Fig. 12. Interestingly, the g factor shows a strong 
dependence on k z ; displaying an oscillatory variation and van- 
ishing intermittently for most subbands. 



top or side gate close to the nanowire or quantum point con- 
tact device. Application of a gate voltage essentially tunes the 
Fermi level (top gate) or the effective lateral confinement (side 
gate) of the nanowire or quantum point contact device such 
that different parts of the energy spectrum fall within the small 
energy window around the Fermi level in which conduction 
is possible. In conventional 77,78 zero-field quantized conduc- 
tance spectra, the conductance changes monotonically in steps 
of 2e 2 /h as a function of gate voltage. This corresponds to 
the case of subband edges at k z = entering the energy win- 
dow which is allowed for transport around the Fermi level. A 
hole nanowire can display quite a different conductance spec- 
trum due to its strong non-parabolic dispersion. For a hole 
nanowire with the dispersion show in Fig. 11(a), the conduc- 
tance can both increase and decrease in steps of 2e 2 /h and 
also 2 x (2e 2 /h). This behavior is caused by off-center sub- 
band maxima, examples of which are seen at points A and B 
in the hole dispersion shown in Fig. 1 1(a), as well as switches 



V. DISCUSSION AND CONCLUDING REMARKS 

We have calculated Zeeman spitting of states in hole 
nanowires, subject to a magnetic field B z pointing parallel 
to the wire axis (the z direction) and defined by three dif- 
ferent types of confinement (cylindrical hard-wall, cylindrical 
soft-harmonic, and rectangular hard-wall). Our main focus 
has been to elucidate the properties of zone-center subband- 
edge states, although a few results are also given for hole- 
wire states having finite wavevector k z for motion along the 
wire. We have disentangled spin-splitting contributions aris- 
ing from (a) the bulk-hole Zeeman effect in the semiconductor 
material, and (b) the coupling of the applied magnetic field to 
the orbital wire bound states. We termed the contribution (a) 
the bulk-material contribution to the g factor and indicated it 
by the symbol g%, where the index a labels the quasi- ID hole- 
wire subbands. The orbital contribution (b) turned out to have 
two separate terms. The one indicated by g" rh diag reflects the 
direct coupling of hole orbital angular momentum component 
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L z with the magnetic field and displays HH-LH splitting , and 
the second one, labelled g" rb mix , arises from HH-LH mixing. 

Quantum confinement of holes in a wire renders their eigen- 
states to be coherent superpositions of HH and LH com- 
ponents. Although these are pure states in the quantum- 
statistical sense, the spin polarization of such HH-LH mixed 
states is not fixed to be 3/2 as a simple analogy with spin- 1/2 
systems would suggest. We have used the recently discussed 63 
invariant-multipole decomposition of the spin-3/2 density ma- 
trix to universally characterize spin polarization of hole-wire 
subband edge states. Within limits, the properties of low in- 
dex (highest-in-energy) hole subbands turned out to be quite 
robust with respect to changes in wire cross-sectional shape 
and confinement strength. Also, many qualitative features are 
already captured within models where band-warping due to 
cubic crystal symmetry is neglected. 

HH-LH mixing is illustrated by radial profiles of the spin- 
3/2 dipole moment (i.e., spin polarization) and turns out 
to have a profound impact on spin splitting. Variations in 
the spin-polarization profiles between different subband-edge 
states are reflected in fluctuations of the individual g-factor 
contributions, as well as the total nanowire-hole g factor. It is 
possible, by addressing different hole-subband edges, to study 
nontrivial hole spin-polarization states that have no counter- 
part in spin- 1/2 systems, such as zero-polarizations states that 
have a large spin-3/2 octupole component. The bulk-material 
contribution to the hole-wire g factor could be investigated 
separately in wires made from diluted-magnetic semiconduc- 
tor materials, where its magnitude can exceed that of the or- 
bital terms by a few orders of magnitude. 

The spin-splitting properties of individual hole-wire sub- 
band edges turn out to be tunable by changing the aspect ratio 
in a rectangular-crossection wire geometry. This is of inter- 
est to applications involving electrostatically defined quantum 
point contacts or gated self-assembled wires, where the wire 



shape can be manipulated in situ. 

Our calculations were performed using Luttinger parame- 
ters for GaAs. However, results obtained within the spherical 
approximation will apply to other materials having a similar 
ratio 7s/7i- The exact quantitative results will probably be 
different once cubic corrections are taken into account. How- 
ever, as our comparison between fully cubic and approximate 
spherical models has shown, qualitative trends for low-index 
hole-wire subband edges are already captured in the spherical 
model. Hence, we believe that our results can be applied to 
materials other than GaAs. 

The effects of linear-in-wavevector spin-orbit couplings has 
not been considered. The relative strength of bulk-inversion- 
asymmetry (BIA) induced k-linear terms, as compared with 
leading-order terms in the Luttinger Hamiltonian, depends 
strongly on the material and the direction of motion with re- 
spect to crystallographic axes. 39 Given that cubic corrections 
have not been too important for low-lying hole-wire subband 
edges, we expect the same to hold for BIA terms. Also, while 
previous studies have shown that, in principle, electric fields 
applied perpendicular to nanowires affect hole g factors 79 (es- 
sentially because the HH-LH mixing changes), this effect be- 
comes relevant only at quite large field magnitudes that are not 
likely to exist in typical point-contact or nanowire samples. 
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